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Abstract. The observed splittings of solar oscillation fre- 
quencies can be employed to separate the effects of inter- 
nal solar rotation and to estimate the contribution from a 
large-scale magnetic field or any latitude-dependent ther- 
mal perturbation inside the Sun. The surface distortion 
estimated from the rotation rate in the solar interior is 
found to be in good agreement with the observed oblate- 
ness at solar surface. After subtracting out the estimated 
contribution from rotation, there is some residual signal in 
the even splitting coefficients, which may be explained by 
a magnetic field of approximately 20 kG strength located 
at a depth of 30000 km below the surface or an equiv- 
alent aspherical thermal perturbation. An upper limit of 
300 kG is derived for a toroidal field near the base of the 
convection zone. 

Key words: Sun: activity - Sun: interior - Sun: magnetic 
field - Sun: Oscillations 



1. Introduction 

Rotational splittings of solar oscillation frequencies have 
been successfully utilized to infer the rotation rate in 
the solar interior. To first order, rotation affects only the 
splitting coefficients which represent odd terms in the az- 
imuthal order m of the global resonant modes. The even 
terms in these splitting coefficients, which reflect the Sun's 
effective acoustic asphericity, can arise from second order 
effects contributed both by the rotation and magnetic field 
as also from latitudinal temperature variations. Since the 
rotation rate can be inferred using the odd splitting coeffi- 
cients, the inferred profile can be used to estimate the sec- 
ond order effects. These can then be subtracted from the 
observed even coefficients to estimate the magnetic field 
strength (Gough & Thompson 1990) or other latitudinal 
variations in sound propagation speed. The distortion in- 
troduced by rotation can be compared with the measured 
oblateness at the solar surface. 

The even coefficients of splittings are fairly small, and 
no definitive results have so far been obtained regarding 
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the magnetic field strength in the solar interior. With the 
good quality data now becoming available from GONG 
(Global Oscillation Network Group) and MDI (Michclson 
Doppler Imager) projects, it is desirable to investigate the 
possibility of inferring the strength of magnetic field in so- 
lar interior. There should also be some shift in the mean 
frequency for each n, t multiplet due to second order ef- 
fects from rotation and magnetic field, which can also be 
estimated. It is difficult to measure this frequency shift 
from observed data as it is hard to separate it from the 
effects of other uncertainties in the spherical structure of 
the Sun. Nevertheless, these frequency shifts can affect 
the helioscismic inferences and it would be interesting to 
estimate their effect. 

2. The technique 

The frequencies of solar oscillations can be expressed in 
terms of the splitting coefficients: 
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where "Pj (m) are orthogonal polynomials of degree j in m 
(Ritzwoller & Lavely 1991; Schou, Christensen-Dalsgaard 
& Thompson 1994). The odd coefficients a\, as, as, . . . can 
be used to infer the rotation rate in the solar interior, 
while the even coefficients arise basically from second or- 
der effects due to rotation and magnetic field. Since forces 
due to rotation or magnetic field in the solar interior are 
smaller by about 5 orders of magnitude as compared to 
gravitational forces, it is possible to apply a perturba- 
tive treatment to calculate their contribution to frequency 
splittings. In this approach, we estimate the effects of ro- 
tation and magnetic field on the frequencies but without 
explicitly constructing a model of a rotating, magnetic 
star. 

We adopt the formulation due to Gough & Thompson 
(1990), with the difference that we include perturbation 
in the gravitational potential and also assume differential 
rotation in the interior, though the symmetry axis of mag- 
netic field is taken to coincide with rotation axis. 
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In an inertial frame the oscillation equations can be 
formally written as 

+ P uj 2 i = cjM£ +Af£ + B£, (2) 

where 

C£ = V(pc 2 V-£ + £-Vp)-(V-£ + £-Vlnp)Vp 
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Here Bi = V x x B) is the linearized Eulerian pertur- 
bation to magnetic field, B, £ is the displacement eigen- 
function, v = f2 x r is the velocity due to rotation, and 
p, p, c s are respectively, pressure, density and sound speed 
in the equilibrium state. 

In the presence of rotation and magnetic field the 
equilibrium state will naturally undergo a distortion that 
needs to be included in the calculations. To account for 
this deformation we consider a transformation to map each 
point r in the distorted star to a point x in the spherical 
volume occupied by the undistorted star by a transforma- 
tion 



x= (1 + ha(r) + h B (r))r, 



where the functions hn(r) and /is(r) which depend on the 
rotation and magnetic field respectively, are to be deter- 
mined by solving the equations for equilibrium in a dis- 
torted star (Gough & Thompson 1990). This will give us 
the perturbation to a nonrotating spherically symmetric 
solar model and the extent of distortion at the surface 
may be compared with observed values. Here, x is chosen 
so that x = R can be regarded as the distorted solar sur- 
face, where R is the radial distance of the outermost layer 
included in the solar model. Similarly, various equilibrium 
quantities are also expressed in the form 



p(r) = p {x) + po(x) + p B (x). 



(8) 



In all these expansions higher order terms have been ne- 
glected. 

We consider the terms on the right hand side of Eq. 0, 
as perturbations to basic equations for linear adiabatic 
oscillations for non-magnetic and non-rotating star. Ro- 
tation introduces a first order perturbation through M. 
which gives the odd splitting coefficients, while magnetic 
field can only give rise to even terms in m and contributes 
to the even splitting coefficients. The distortion from a 
spherically symmetric equilibrium state also introduces 
even order terms. The relative magnitude of contributions 
from rotation and magnetic field will, of course, depend 



on the rotation rate and magnetic field strength. For the 
solar case we know that odd splitting coefficients arising 
from the first order effect of rotation are much larger than 
the even coefficients and we therefore expect the mag- 
netic field to make a comparatively smaller contribution. 
We must therefore include the effect of rotation to second 
order, while magnetic field and distortion effects need be 
retained only to first non- vanishing terms. The first order 
perturbation arising in frequencies on account of rotation 
also introduces a perturbation to eigenfunctions which will 
give a second order contribution. We can formally express 
the frequency and eigenfunction as 

u> = loq + GJi + uj 2 , £ = £o+£i- ( 9 ) 
Retaining terms to second order, we get 

£o(£ + €i) + £n£o + £b£o + poH + 2^i)(£ + Si) 

+p (ujf + 2w ^2)£o + Pn^o + pbojqZo 
= w M{t + t 1 )+wiMt +Arta+Bt . (10) 

Here, £r> and Cb are the perturbations to C arising from 
distortion of equilibrium state due to rotation and mag- 
netic field respectively. Taking the scalar product with £q 
and integrating over the entire volume, we recover 

+(Co(l3 - C B - PbuDZq) - ^(Mi^o) 
-2woWi<po€S{i) +wi(£S^£o) +^{C a M^), (11) 



(7) where the angular brackets denote 



(f(x, 0, </>)) = / f(x, 0, (j))x 2 sin dx d9 d<j) 
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The first order correction to frequency is given by 



2<M^o> ' 



(12) 



(13) 



while perturbation to the eigenfunction may be calculated 
using 



£f i + po^oii = -2p woWi£o + uoM£ 



(14) 



The observed odd splitting coefficients can be used 
to infer the rotation rate inside the Sun (Thompson et 
al. 1996; Schou et al. 1998). We approximate this rotation 
rate using the first three terms in the expansion of the 
angular velocity, 



n{r, 0) = fio(r) + n 2 (r) cos 2 + Q 4 (r) cos 4 0, 



(15) 



where is the colatitude. This rotation rate is then used 
to compute the second order rotational contribution to 
frequency splitting, which may be subtracted from the ob- 
served splittings to obtain the residual which may be due 
to magnetic field, any other velocity field or asphericity in 
solar structure. 
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In the present analysis we use only the toroidal mag- 
netic field, taken to be of the form, 
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n n t \ dPk ( cos ' 



(16) 



with the axis of symmetry coinciding with the rotation 
axis. Here Pk{x) is the Legendre polynomial of degree k. 
The Lorentz force due to a field of this form can be written 
as 
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Each of this term can be treated separately and the results 
can be combined to yield the net effect. 

We calculate the second order frequency shift due to 
rotation and magnetic field for each value of m and then 
use Eq. [l]to obtain the corresponding splitting coefficients. 
These can then be compared with observed coefficients 
from GONG (Hill et al. 1996) or MDI (Rhodes et al. 1997) 
data. To evaluate the angular integrals we use the follow- 
ing recursion relations 



cos 6YJ 
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Since we have used only the first two terms in the expan- 
sion of rotation rate as a function of latitude, we restrict 
to calculation of the splitting coefficients 02 and 04 in this 
work. 

3. Results 

We use the rotation rate inferred from the GONG data for 
the months 4-14 (Antia, Basu & Chitre 1998) to estimate 
the second order frequency shift and the corresponding 
splitting coefficients 02 and 04, as outlined in the previous 
section. We incorporate all the second-order contributions 
arising from rotation, including those from the distortion 
of equilibrium state and the perturbation to the eigen- 
functions. Although there may be some variation in rota- 
tion rate with time, the estimated variation is very small 
and its effect on the inferred splitting coefficient would be 
much smaller than the errors in observed values. 



3.1. Shift in the mean frequency 

In principle, the shift in the mean frequency arising from 
second order effects of rotation can be calculated with 
the help of the prescription outlined in the previous sec- 
tion, by taking the spherically symmetric component of 
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Fig. 1. The shift in mean multiplet frequencies due to 
second order effects from rotation is shown by the crosses, 
while the squares (points with positive 8v) show the fre- 
quency shift due to general relativistic effects. 



the perturbing force (A = term in Eq. [17]). However, 
this will also change the mass, radius and luminosity of 
the solar model. The change may be smaller than the er- 
rors in observed radius or luminosity, but it may tend to 
give a different estimate for modified frequency compared 
to what will be obtained if the observed constraints on 
mass, radius and luminosity were to be exactly applied. 
Hence, for obtaining a consistent estimate of the effect 
of distortion, we construct a spherically symmetric solar 
model with correct mass, radius and luminosity by modi- 
fying the effective acceleration due to gravity, g to account 
for the spherically symmetric component of forces due to 
rotation. The difference in frequency of this model in rela- 
tion to a standard, non-rotating model would give the fre- 
quency shift due to distortion. All the other second-order 
rotational terms are added to this shift, to obtain the to- 
tal shift in frequency due to rotation which is displayed in 
Fig. 1. This figure includes all modes with 0.5 < v < 4.5 
mHz and £ < 250. The corrections to mean frequencies 
due to general relativistic effects as discussed towards the 
end of this subsection, are also shown in the figure. 

This relative frequency shift, which is less than 10~ 5 , is 
nonetheless comparable to the estimated errors in the ob- 
served frequencies and the correction should, in principle, 
be applied while doing inversions (e.g. Gough et al. 1996) 
for the Sun's spherical structure. In order to estimate the 
error introduced by neglect of this effect, we can carry out 
an inversion for sound speed and density in the solar inte- 
rior using this frequency shift due to rotation as the fre- 
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Fig. 2. The correction to sound speed and density as 
inferred from hclioseismic inversions, arising from the fre- 
quency shifts shown in Fig. 1. The thick continuous and 
short-dashed lines show, respectively, Sc 2 /c 2 and Sp/p due 
to frequency shifts from rotation (crosses in Fig. 1), while 
the thin lines show the same arising from combined fre- 
quency shifts due to rotation and general relativity. 



quency difference and the results are shown in Fig. 2. The 
inversions are performed using a regularized least squares 
inversion technique (Antia 1996). The resulting 6c 2 /c 2 and 
Sp/p are almost an order of magnitude less than the esti- 
mated errors in inversions. 

As an aside, we note that the internal rotation rate 
from Antia et al. (1998) adopted in our study was ob- 
tained assuming a spherically symmetric background state 
for the Sun, as is usual for inversions for the solar rotation. 
We realise that both the mean frequencies of solar oscil- 
lations and the rotational splittings will be modified by 
departures in the equilibrium solar model from spherical 
symmetry, as discussed in this paper. In order to estimate 
the resulting shift in rotational splittings we would need 
to calculate the third order terms in perturbation expan- 
sion of Gough & Thompson (1990). We have not included 
these terms in our analysis, but we expect that their con- 
tribution would have the same relative magnitude of 10 
as that found for the shift in mean frequencies. This is 
clearly, much smaller than the estimated errors in split- 
ting coefficients in current helioseismic data sets. There- 
fore we do not expect the rotational splittings and hence 
the inverted rotation rate to be significantly affected by 
this higher order effect. 

It may be noted that mean frequencies of f-modes get 
diminished by up to 15 nHz on account of the effect of 
rotation. Since rotation effectively reduces the acceleration 
due to gravity g, this leads to a decrease in the frequencies 
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Fig. 3. The relative shift in mean multiplct f-mode fre- 
quencies due to rotation 

of f-modes. The relative change in f-mode frequencies is 
shown in Fig. 3. If this effect is taken into account the 
estimated solar radius using f-mode frequencies (Schou et 
al. 1997, Antia 1998) would effectively be decreased by 
about 4 km. This is again much less than the systematic 
errors in estimated radius, though the decrease is larger 
than the statistical errors (Tripathy & Antia 1999). 

It is interesting to note that apart from second or- 
der effects of rotation, there would also be corrections to 
the frequencies arising from general relativity. The rela- 
tivistic effect can be measured by Gm(r)/ (rc 2 ), where G 
is the gravitational constant, m(r) is the mass contained 
within spherical shell of radius r, and c the speed of light. 
Fig. 4 shows this ratio in a solar model and it can be seen 
that it is comparable to the ratio of centrifugal to grav- 
itational forces. It is possible to calculate a solar model 
using Oppcnhcimcr-Volkoff equation of relativistic stellar 
structure instead of the standard equation of hydrostatic 
equilibrium: 

dp _ G(p + p/c 2 )(m + 4:Trr 2 p/c 2 ) 
~cfr ~ r 2 (l -2Gm/rc 2 ) ' 



(21) 



It is clear that general relativistic effect would be of op- 
posite sign to that due to rotation, as rotation effectively 
reduces the acceleration due to gravity, g, while the rel- 
ativistic correction tends to increase it. Thus there is a 
partial cancelation between the two effects. It is possible 
to calculate the change in solar models due to the rela- 
tivistic effect, although a detailed calculation of frequen- 
cies using relativistic stellar oscillations equations would 
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Fig. 4. The ratio of centrifugal force and the gravitational 
force (continuous line) in a solar model is shown as a func- 
tion of the radial distance r. The dashed line shows the 
ratio Gm(r)/(rc 2 ) which gives the effect due to general 
relativity. 



to those expected from rotational distortion. There may 
be some residual arising from other effects, like magnetic 
field or other asphcricitics. The contribution from mag- 
netic field is indeed expected to vary with solar cycle and 
may account for the variation in P4 component, if the 
variation is in fact real. 

It is also possible to estimate the global parameters for 
the Sun, like angular momentum, rotational kinetic energy 
and gravitational quadrupole and hexadecapole moments 
due to rotational distortion (Pijpers 1998) and the results 
are summarized below: 

Moment of Inertia, / = 7.11 x 10 53 gm cm 2 , (22) 
Angular Momentum, H = 1.91 x 10 48 gm cm 2 s~\ (23) 
Kinetic Energy, T = 2.57 x 10 42 gm cm 2 s~ 2 , (24) 
Quadrupole Moment, J 2 = -2.18 x 1(T 7 , (25) 
Hexadecapole Moment, J 4 = 4.64 x 1(T 9 , (26) 

which are consistent with estimates of Pijpers (1998), who 
obtained his estimates by working in terms of kernels for 
the various quantities. The value of J 2 will yield a pre- 
cession of the perihelion of planet Mercury by about 0.03 
arcsec/century, which is small enough to maintain consis- 
tency of the general theory of relativity. 



require considerable effort and is beyond the scope of the 
present work. To a first approximation we may calculate 
the effect by using the normal equations of stellar oscil- 



lations with gravity modified according to Eq. 21. Such a 



calculation shows that the effect of relativity more or less 
cancels the frequency shift due to rotation for low degree 
modes. The frequency shift due to general relativity are 
also shown in Fig. 1. If this frequency shift is added to 
the contribution arising from rotation then the effect on 
helioseismic inversion is significantly reduced in the solar 
core as can be seen from Fig. 2 (compare the thick and 
thin lines). 

3.2. Oblateness due to rotation 

During the course of computing the splitting coefficients, 
it is necessary to calculate the deformation induced by ro- 
tation as outlined by Gough and Thompson (1990). This 
deformation may be compared with the observed oblate- 
ness at the solar surface. The surface amplitudes of the 
P2{cos8) and P 4 (cos6') components of deformation are 
found to be — 5.84x 10 -6 and 6.2x 10~ 7 respectively, which 
are consistent with the estimates obtained by Armstrong 
& Kuhn (1999). These can be compared with measured 
values of -(5.44 ± 0.46) x 10~ 6 and (1.48 ± 0.58) x 10" 6 
respectively, from MDI measurement during 1997 (Kuhn 
et al. 1998). Kuhn et al. (1998) find a large temporal vari- 
ation in the P4 component, but it is not clear if the vari- 
ation is statistically significant. It can be seen that the 
measured values of solar oblateness are reasonably close 



3.3. Second order splitting due to rotation 

The contribution to splitting coefficients a-i and 04 due 
to rotation is shown in Fig. 5. This contribution needs 
to be subtracted from the observed splitting coefficients 
for obtaining the residual contribution which may arise 
from effects due to magnetic field, other velocity fields 
or asphericity in solar structure. Since the errors in in- 
dividual splitting coefficients are too large to give signifi- 
cant differences, we average over 30 neighbouring modes in 
w = vj (£+1/2) and the corresponding results are shown in 
Fig. 6. There is reasonable agreement between the GONG 
data for months 4-14 (23 August 1995 to 21 September 
1996) and MDI data for the first 360 days of its opera- 
tion (1 May 1996 to 25 April 1997). It is well known that 
the even splitting coefficients vary with solar activity cy- 
cle (Libbrecht & Woodard 1990; Dziembowski et al. 1998; 
Howe, Komm & Hill 1999) and there may not be agree- 
ment between observations taken at different epochs. But 
in the present case there is considerable overlap in period 
and the observations are near the minimum phase of solar 
activity, when these coefficients are not expected to vary 
significantly. 

The difference between the observed splitting coeffi- 
cients and the estimated contribution from rotation is sig- 
nificant for modes with turning points in the convection 
zone. For modes penetrating more deeply, the errors are 
larger and the difference is probably not significant. 
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Fig. 5. The splitting coefficients and 0,4 from effects of rotation. 



40 



20 



X 



20 



□ GONG 
x VfDl 

■ Rotation 



{ f T i ■ lis 1 




t 1 . itilM J 



X 



0.2 



0.4 



0.6 



0.8 



r./R„ 




Fig. 6. The splitting coefficients 02 and 04 plotted as a function of the lower turning point for the modes. The 
coefficients from GONG and MDI data are compared with the contribution from rotation. Each point represents an 
average over 30 neighbouring modes. 



3.4. Splitting due to magnetic field near the base of the are expected from such a held by assuming the magnetic 
convection zone field to be given by Eq. |l| with 



l(r) = ( VS^CRl - (^) 2 ) if |r - r | < d (27) 
[ otherwise 



There have been some suggestions that a significant 

toroidal magnetic field may be concentrated in a layer where po is the gas pressure, /3o is a constant giving the 

around the base of the convection zone (Dzicmbowski & ratio of magnetic to gas pressure, r and d are constants 

Goode 1992). We therefore first investigate splittings that defining the mean position and thickness of layer where 
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the field is concentrated. Fig. 7 shows the splitting co- 
efficients resulting from a toroidal magnetic field of this 
form concentrated at the base of the convection zone 
(r = O.713i? and d = 0.02R Q ). The splitting shown in 
this and subsequent figures includes both the direct and 
distortion contributions as defined by Gough and Thomp- 
son (1990). 

The coefficients a 2 and 04 from a toroidal magnetic 
field concentrated near the base of the convection zone 
have a characteristic signature for modes with turning 
point near the base of the convection zone; it should be 
possible to detect such a signal in the observed splittings 
if a strong enough magnetic field is indeed present in 
these layers. The computed splittings, particularly for the 
deeply penetrating modes in Fig. 7 show a great spread, 
which is characteristic of the splittings arising from a thin 
magnetic layer. We return below to the use that can po- 
tentially be made of this signature. In the present study, 
however, we choose to average over neighbouring modes, 
as discussed in Section 3.3, which suppresses this spread. 
Our rationale is that the errors in the real data are too 
large for the spread to be visibly distinguished from noise 
in the measured splittings at present. Thus we take aver- 
ages over neighbouring modes and compare the residual 
after removing the contribution due to rotation with the 
expected splitting from the magnetic field and the results 
are shown in Fig. 8. Note that even after averaging a clear 
signature of the magnetic field is seen in the splitting co- 
efficients. Since we are comparing the average over the 
same set of modes for the observed splittings and com- 
puted splittings for magnetic field, we should be able to 
get some estimate of magnetic field if a strong enough field 
does indeed exist. From Fig. 8 it can be seen that there 
is no clear signature of any feature near the base of the 
convection zone in the observed splittings, and hence we 
can only set an upper limit on the magnetic field in this 
layer. This will of course, depend on the thickness of the 
magnetic layer. Since there is no clear signature of any 
signal near the base of the convection zone, for quantita- 
tive purpose we take the difference between the lowest and 
highest point in the range 0.6 < r t /i? Q < 0.8 in observed 
splitting coefficients. For a 2 this difference is 8.7 nHz for 
MDI and 7.0 nHz for GONG data, while computed split- 
tings with /3o = 10~ 4 show a difference of 12.6 nHz for a 
half-thickness of 0.02i? Q . Thus, we can put an upper limit 
of 0.7 x 10~ 4 on /3o which corresponds to a magnetic field 
strength of 300 kG for a layer of half-thickness 0.02i? Q 
near the base of the convection zone. Similar analysis for 
splitting coefficient 04 yields a slightly larger upper limit 
of 400 kG. These limiting values are close to what was 
obtained by Basu (1997) using a similar technique and is 
also consistent with the value independently inferred by 
D'Silva & Choudhuri (1993). Note, this limit roughly in- 
creases as 1/ y/d, and clearly, if the thickness of this region 
is smaller, the upper limit would be larger. It should be 
noted that the tachocline, where the rotation rate under- 



goes a transition from differential rotation in the convec- 
tion zone to a solid-body like rotation in the radiative inte- 
rior may have a thickness as small as 0.01i? Q (Basu 1997; 
Antia, Basu & Chitre 1998). With this thickness the upper 
limit on magnetic field would naturally be increased. 

There is the possibility of distinguishing seismologi- 
cally between magnetic layers of different thicknesses by 
using modes that penetrate well beneath the magnetic 
layer. A thin layer will induce a signature in the ai and a 4 
coefficients which is periodic in mode frequency (Gough 
& Thompson 1988; Vorontsov 1988; Thompson 1988), 
in much the same way that the rather sharp transition 
near the base of the convective envelope produces a peri- 
odic signature in the mean frequencies (e.g., Gough 1990; 
Basu, Antia & Narasimha 1994; Monteiro, Christensen- 
Dalsgaard & Thompson 1994). Indeed it is this signa- 
ture which is largely responsible for the vertical spread of 
points for modes with turning points at radii r ^ 0.6i? in 
Fig. 7. Basu (1997) attempted to use this oscillatory signal 
to obtain an upper limit on magnetic field near the base 
of convection zone (see also Gough & Thompson 1988). 
Fig. 9 shows la<i for modes with I < 10 for magnetic field 
concentrated near the base of the convection zone, with 
two different values of d. It is clear that the amplitude 
of oscillatory signal varies significantly with d. However, 
the observed splitting coefficients for low values of i have 
large errors and it is difficult to extract the small oscilla- 
tory signal from these. 

3.5. Field in the upper convection zone 

Having considered a magnetic field at the base of the 
convection zone, where theory suggests a field might be 
stored, we consider where else the data might indicate 
the presence of magnetic field. There is no signature for 
the presence of significant magnetic field in the radiative 
interior, since the averaged residual splitting after correct- 
ing for rotation seem to be consistent with zero. However, 
within the convection zone there is some significant resid- 
ual splitting, which could be due to the effect of a magnetic 
field. An inspection of these residuals indicates the exis- 
tence of a peak around r — 0.96i? Q , and indeed, if it is 
due solely to magnetic field, the field may be distributed 
around this depth (w 28000 km). It may be noted that 
this is approximately the depth to which shear layer seen 
in rotation profile extends (Antia, Basu & Chitre 1998; 
Schou et al. 1998). 

We now attempt to estimate splittings due to the field 
concentrated in this region. Fig. 10 shows the splittings 
due to a few magnetic field configurations which are con- 
centrated in the upper part of the convection zone. A com- 
parison of these with the observed splittings indicates that 
there may be an azimuthal magnetic field with (3 < 10 -4 
(i.e., B w 20000 G), with peak around r = 0.96i? Q . 

The possible existence of a magnetized layer with field 
of order 20 kG located around r — 0.96i? Q is, indeed, a 
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Fig. 7. The splitting coefficients a 2 and 04 from a toroidal magnetic field concentrated near the base of the convection 
zone, plotted as a function of the lower turning point for the mode. Magnetic field is given by Eqs. (16,27) with k = 2, 
/?o = 1CP 4 , r = O.713i? (shown by the vertical line in the figure) and d — 0.02i? Q . 
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Fig. 8. The splitting coefficients 02 and 04 from a toroidal magnetic field concentrated near the base of the convection 
zone, plotted as a function of the lower turning point of the modes. Each point represents an average over 30 neigh- 
bouring modes. The estimated contribution from rotation has been subtracted from the observed splittings plotted in 
the figure. Magnetic field is given by Eqs. (16,27) with k = 2, /3 = 1(T 4 , r = 0.713^0 and d = 0.02.R©. 



significant inference drawn from our analysis. The physical 
interpretation for the origin of such a moderately strong 
magnetic field at this depth below the Sun's surface is 
naturally a challenging task for theories of solar dynamo 
to accommodate. It may be useful to recall here that 



the numerical simulations of the Sun's outer convection 
zone (Nordlund 1999) indicate a major presence of down- 
ward moving plumes. It is conceivable that these down- 
drafts could gather the turbulent magnetic field in the sub- 
surface layers and carry them to depths in the convective 
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Fig. 9. The splitting coefficient ct2 from a toroidal magnetic field concentrated near the base of the convection zone, 
plotted as a function of the frequency for modes with I < 10. Magnetic field is given by Eqs. (16,27) with k = 2, 
/3q = 1CP 4 , ro = O.713i?0 and the value of d as marked in the panels. 
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average over all f-modes. The estimated contribution from rotation has been subtracted from the observed splittings 
plotted in the figure. Magnetic field is given by Eq. with /3 = 10~ 4 , and the value of r , d and k as marked in the 
figure. 



envelope until some sort of equipartition is reached. Inter- 
estingly, the density, p at a depth of 25-30 Mm is upwards 
of 4 x 10~ 3 gm cm' 3 , while the downward velocity for the 
plumes is of order 500 m s _1 . The dynamical pressure of 
the plumes, pv 2 ^ 10 7 dyne cm~ 2 , then becomes compara- 



ble with the magnetic pressure, B 2 /87r, corresponding to 
a field strength of 20-30 kG. It is, therefore, tempting to 
envisage the formation of such a magnetized layer by the 
pounding of the downdrafts which tend to concentrate the 
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field at depths where the equipartition of the kind outlined 
above is approached. 

In this study we have assumed a smooth toroidal mag- 
netic field, but in practice we do not expect such a field 
inside the convection zone. Turbulence may be expected 
to randomize the magnetic field and such a field may not 
be expected to produce any significant distortion in the 
equilibrium state. The direct effect of magnetic field will 
still be felt though the contribution would be different. 
Thus our results may be treated as indicating the order of 
magnitude of field that may be expected if the observed 
splitting coefficients are indeed due to the magnetic field. 
If the field is concentrated in flux tubes which occupy only 
a small fraction of the volume then the required magnetic 
field could be correspondingly larger. If we assume that 
the flux tubes occupy a fraction / of the total volume, 
the magnetic field strength should increase by 1/V7- If 
we consider only direct contribution to the splittings then 
it turns out that a-i is always negative for all toroidal field 
configurations that we tried and hence such a contribu- 
tion is not likely to explain the observed splittings. But 
a different magnetic field configuration, e.g., poloidal field 
might produce a 2 with the required sign using only di- 
rect contribution. The order-of-magnitude splitting caused 
by a magnetic field is la-ijv ~ /3q ~ v 2 A /c 2 s , where va is 
the Alfven speed. We therefore regard it as unlikely that 
a different magnetic field configuration would produce a 
markedly different answer for the field strength required 
to account for the observed signal in ai and a^. 

A nonmagnetic latitudinally-dependcnt perturbation 
to the wave propagation speed might be responsible for 
the signal we detect (cf., Gough & Zweibel 1995). Once 
again we may expect a perturbation of order 10 -4 located 
in the region around r — 0.96i?© to yield the observed 
splittings. Gough ct al. (1996) inferred a perturbation of 
that magnitude, of unspecified origin, from earlier GONG 
data. A temperature variation of order 10K, suitably con- 
fined, might conceivably produce a similar signature. In 
fact, Kuhn (e.g., Kuhn 1996) has argued that the ther- 
mal shadow of belts of magnetic flux near the bottom of 
the convection zone can have a significant effect on the 
even a-coefficients. But Kuhn's models show the largest 
temperature perturbations occurring in the very superfi- 
cial superadiabatic layer, at a depth of a small fraction of 
one per cent of the solar radius. Such a perturbation alone 
would be consistent with the f-modes having a small resid- 
ual splitting, but would not explain the apparent overturn- 
ing of the p-modc splittings at r t ~ 0.96i?. A magnetic 
field at some depth below the surface may explain both 
aspects. We certainly do not rule out the possibility that 
some nonmagnetic asphericity, which we have not consid- 
ered in detail in this study, may account for some of the 
observed splittings. 



4. Conclusions 

Second order correction to mean frequencies due to rota- 
tion is comparable to the error estimates in the observed 
frequencies. The error in helioseismic inversion introduced 
by the frequency shift due to rotation is ^ 1CP 4 , which is 
much smaller than the estimated errors in inversions. Fur- 
ther, a part of this frequency shift is expected to be nulli- 
fied by the general relativistic effects. The shift in f-mode 
frequencies due to rotation can reduce the estimated solar 
radius by 4 km. The distortion due to rotation can yield 
surface oblateness of —5.8 x 1CP 6 and 6.2 x 1CP 7 in the 
P2(cos6>) and P4(cos8) components, respectively. This is 
in reasonable agreement with observed oblateness at the 
solar surface (Kuhn et al. 1998) and it appears that most 
of the observed distortion is accounted by the seismically 
inferred rotation rate in solar interior. The quadrupole mo- 
ment J2 = —2.18 x 1CP 7 resulting from rotational distor- 
tion is small enough to maintain consistency of the general 
theory of relativity. 

After subtracting the estimated contribution from ro- 
tation to the splitting coefficients ai and 04 from the ob- 
served splittings, there is a small residual which is statis- 
tically significant in the convection zone. This could arise 
from a magnetic field. From the magnitude of residual in 
observed splittings we can tentatively conclude that mag- 
netic field with P w 10 -4 may be present in the upper 
part of the convection zone. This corresponds to an az- 
imuthal magnetic field of ss 20 kG around r = 0.96i? Q . 
However, we cannot rule out the possibility that this sig- 
nal in splitting coefficients may arise from some aspheri- 
cal perturbation to the temperature field. This would be 
practically indistinguishable from the effect of a magnetic 
field using just the mode frequencies (Gough & Zweibel 
1995); but complementary analyses such as time-distance 
hclioseismology might be able to distinguish them, since 
the local direct effect of a magnetic field on the waves is 
anisotropic, whereas that of a temperature perturbation 
is not. 

A toroidal magnetic field that is concentrated near the 
base of the convection zone gives a characteristic pattern 
in the splittings for modes with lower turning point in that 
region. Since no such signal is seen in observed frequencies, 
we can put an upper limit of about 300 kG on the strength 
of the magnetic field in this region. 
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